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Gas Phase Recombination of Hydrogen and Deuterium Atoms 

Daniel W. Tralnor* , David 0. Ham**, and Frederick KaufYnan 
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Pittsburgh* Pennsylvania 15213 


ABSTRACT 

Rate constants for the reaction H + H + M ■+ H 2 + M, with M = H 2 , 

He, and Ar were measured over the temperature range 77 to 298°K. Hydrogen 
atoms were produced by thermal dissociation and absolute atom concentrations 
were measured through use of self-balancing, isothermal catalytic probe 
detector. The specific rate constants were 8.1 + 0.4 x 10” 33 , 7.0 ±0.4 x 10“ 33 , 
and 9-2 + 0.6 x 10 ” 3 3 cm 6 molecules -2 sec -1 at 298°K for M = H 2 , He, and Ar 
respectively; these values rising to 18. 5 ± 2.2 x 10 ” 33 , 12.0 + 1.5 x 10 " 33 , 
and 27.4 4 , 4.6 x 10 “ 33 cm 6 molecules ” 2 sec -1 at 77°K. For the equivalent 
deuterium atom process with D 2 as the third body, the rate constants are 
6.1 ± 0.3 x 10 “ 33 -ton 6 molecules ” 2 sec ” 1 at 298°K and 15. 1 ± 1.0 x 10 “ 33 cm 6 
molecules " 2 sec ” 1 at 77°K. ,; 3hese values ■ are' conpared with previous experi- 
mental measurements and with recent 5 theoretical calculations. 
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I. Introduction 

As the simplest termolecular reaction and one that is well suited 
for rigorous calculations, the gas phase recombination of hydrogen atoms in 
the presence of third body atoms or molecules, M, continues to generate 
wide interest both theoretically and experimentally. Although there have 
been many experimental determinations of the hydrogen atom recombination 
rate constant 1-16 ranging from the earliest investigation by Smallwood 1 in 
1929, to recent studies of Larkin and Thrush, 10-11 and of Bennett and 
Blackmore, 12-14 its value has remained surprisingly uncertain, there is 
very little known about its temperature dependence , particularly below 
room temperature, and about the effectiveness of different M. Recent room 
temperature studies for M = H 2 , since i960, show poor agreement as they 

report specific rate constants ( d ^ 2 - / [H] 2 [M]) covering a range from 
3-6 to 200 x 10“ 3 3 cm 6 molecule -2 sec -1 . A critical review of these data, 
as well as those obtained in flames and shock tubes, can be found in a 
new .book by Baulch, Drysdale, Lloyd, and' Home. 17 

Aside from its continuing challenge to the experimentalist and the 
extreme dearth of measurements of the temperature dependence and of the M 
effect, the problem has become much more interesting due to the recent cal- 
culations of a quantum mechanical orbiting resonance theory by Roberts, 
Bernstein and Curtiss, 18 and a classical phase space theory by Shul , 
Appleton and Keck. 1 9 > 20 These theories predict quite different behavior 
in the low temperature region, and thus have stimulated much interest in 
experimental studies covering a wide temperature range for a variety of 
third bodies. 

In this paper, we report rate constants for the reaction measured 
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in a flow system over the temperature range 77 to 298°K, and for M = H 2 , 

He, and Ar. As the reaction is second order in H, the absolute H-atom 
concentration must be measured accurately. In addition, experimental con- 
ditions should be selected so as to provide a manageable kinetic model 
for treating the data. For these measurements, an isothermal catalytic 
probe atom detector was utilized as the absolute atom measuring device. 

Its accuracy as an atom detector is dependent, however, both on complete 
atom recombination occuring on the probe and on quantitative energy 
accommodation of the heat released on recombination. We have measured the 
efficiency of our probe for both these effects and found it to be complete. 
The present data thus represent the first extensive and self-consistent 
experimental study of the temperature, isotope, and M-effect at and below 
300°K and can form the basis for fruitful comparison with theory. 

II. Experimental 

A . Flow System 

The apparatus is shown schematically in Fig. 1. Gaseous reactants 
are admitted through needle valves upstream of the reaction tube and their 
flowrates are measured by float type flowmeters (Fisher and Porter) cali- 
brated with a volumetric displacement device (Vol-u-meter, G. K. Porter Co.) 
Calibrations were reproducible to ±. 2 % over a period of two years. The 
reaction tube pressure is measured with a pressure transducer (Validyne 
Engineering Co., Model DP7, ±. 1 psi and Model CD12 Indicator) calibrated 
with a McLeod gauge (Consolidated Vacuum, Model GM-100A) and oil manometer. 

The central 100cm of the quartz reaction tube (1.000 inch i.d., 
total length 180 cm) is surrounded by a copper cooling jacket which, 
in turn, is surrounded by a cylindrical heating mantle, consisting of a long 
central section and two separate end sections, controlled by a solid state 
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temperature controller (HI controls, thermal Model MPKY) . The accessible 
tenperature range is 77 to about 900°K, but the present work covers only 
the lower range, 77 to 298°K. These lower temperatures are obtained by 
passing cold N 2 through the; copper jacket., the gas tenperature being 
measured by three external thermocouples ■( lron-constantan ) placed between 
the reaction tube and the cooling jacket, and by the probe itself cali- 
brated and used as a resistance thermometer. For experiments at 77°K, 
liquid N 2 Is blown through the cooling jacket from a large (160 lit) 
container and collected in small dewar flasks at the exit ports. For 
experiments at temperatures between 77 and 296°K, liquid N 2 is pumped 
through the jacket by suction using a small rotary air purrp (Eberbach Corp.) 
at controlled flow rates. 

B. H-Atom Source 

Hydrogen atoms are produced by the thermal dissociation of H 2 
( 99-999 % pure) on a hot tungsten filament which is suspended in a water 
cooled quartz housing as shown in Fig. 2. H 2 is passed through a liquid 
nitrogen cooled zeolite trap (molecular sieve, 5^, Davison Chemical), through 
the upstream dissoclator, then through a liquid nitrogen cooled glass wool 
trap, before passing over the second dissoclator and through a pyrex to 
quartz seal 14 cm long and 1.2 cm i.d. into the main flow tube. The tungsten 
filaments were prepared by tightly winding 0.015 inch wire on a mandril 
and then stretching it into a loosely wound helix about 30-40 cm long. 

This wire was spot-welded to lead-in rods which are part of a demountable 
O-ring glass joint. This makes it possible to move either dissoclator from 
one position to the other. Initially, a filament was placed in the up- 
stream 'position where it was heated to approximately 2200°K in a continuous 
flow of H 2 for a period of several hours. After this treatment, it is ready 
for use as the atom source in the downstream position and another W-wire 
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assembly is placed in the upstream position. This assures the instant 
availability of a properly pretreated W-filament, and also removes from 
the gas stream any reducible impurity by reaction with H or H2 at high 
temperatures followed by condensation in the packed trap at 77°K. 

C. H-Atom Detector 

The isothermal probe atom detector was constructed of 63 cm of 
Pt/10^Rh thermocouple wire (0.010 inch diameter, Engelhard Industries) 
which was wound as a spiral on a quartz probe support as shown in Fig. 3. 

The probe wire forms one arm of self-balancing Kelvin bridge whose circuit 
diagram is shown in. Fig. 4. A Kelvin rather than Wheatstone bridge network 
is used in order to eliminate or greatly reduce the effect of lead and con- 
tact resistance. This is particularly advantageous under present experimental 
conditions where variable lengths of the copper lead wires are at low 
temperatures and would therefore contribute a changing resistance due to 
their temperature coefficients. 

In normal operation of the self-balancing bridge, the standard 
resistor (0 to 1 ohm, 0.1 ohm steps) is set at> a resistance, Rg, which equals 
that of the probe, Rp, at a sufficiently elevated temperature, AT $ 100°C, 
such that the maximum temperature rise due to atom recombination is less 
than (but not too much less than) AT. Ihe bridge unbalance sigpal, Rg-Rp, 
is amplified by the operational amplifier (Fhilbrick, P85A, gain - 2000) 
and controls the current passed by the transistor (NPN silicon power transistor, 
2N2230). In this manner, since the probe resistance, R D , increases with 
increasing current (temperature) whereas Rg is independent of it, a feedback 
loop is established which keeps the magnitude of Rg-Rp at less than about 
0.001 ohm, corresponding to an uncertainty of less than 0.25°C in the tent- 
perature of the probe wire. 



Ihe probe is completely enclosed in the flow system (see Pig. 1), 
the wires to the Kelvin bridge passing through a vacuum feedthrough. 

Attached to the end of the probe support rod is a small steel cylinder, 
which is surrounded by a solenoid. When the solenoid is activated, the 
probe is thereby magnetically coupled to the external drive mechanism, 
causing the probe to move up-or downstream in the flow tube over a dis- 
tance of 1 meter at any one of ten steady speeds from 0.05 to 2.3 cm se c- 1 
selected by appropriate speed reduction of a synchronous motor (Apcor speed 
reducer. Model 2405 coupled to a Bodine Motor, 1/15 HP, 690 R.P.M., 115 V 
ac 60 CPS). This arrangement of the probe lead wires inside a narrow, 
evacuated tube eliminates the need for sliding rods and 0-ring seals. 

Ihe main gas flow passes througi the punping port at the end of the reaction 
zone, and the effect of the additional gas flow path on conditions in the 
reaction zone is negligible. 

Cylinder gases of the following indicated purity were used without 
further purification: H2 (ultrahigh purity, 99 - 999 %), D 2 (C.P., 99 - 5 %), 

He (high purity, 99 - 995 %), Ar (prepurifled, 99 - 998%) . NO (technical, 98 . 5%) 
and NO Cl ( 97 %) were condensed, pumped and distilled several times before use. 

III. Calibrations and Procedures. 

Catalytic probe atom detectors 21 have been critized as unreliable 
for accurate kinetic work because of their lack of specificity, possible 
inefficiency at collecting and/or detecting atoms, and pertubation of the 
system by acting as an atom sink and thereby producing concentration gradients 
near the probe. Moreover, the general limitations of flow system parameters 
must be kept in mind and checked experimentally . Gaseous impurities were 
minimized through use of a helium leak tight flow system and ultra pure 
.hydrogen. It was further necessary to test whether the probe 
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(a) intercepted and recombined the H-atoms quantitatively; and 

(b) received their total reconciliation energy, i.e. exhibited unit energy 
accommodation. 

A. Atom Collection Efficiency 

To measure the efficiency of the probe as an atom collector, several 
experiments were performed using the set-up shown in Fig. 5- With the probe 
inserted upstream of the H-atom entrance port and placed at position A in 
the main flow tube, addition of small amounts of NO to the atom stream 
permitted measurement of the resulting HNO* chemiluminescence emission 
intensity, I, by the photocurrent of a photomultiplier tube (RCA 7265, S-20 
response) viewing the flow tube through an interference filter (IR Industries) 
which had a transmission maxinium at 6925 ft and a bandwidth (fwhm) of looft. 

Since I = Io[H][NO], 22 the intensity measures relative H-atom concentrations 
when [NO] is held constant. The probe was then moved downstream past the 
H-atom source inlet to position B and the HNO* intensity was re-determined. 
Experiments indicated that the fraction of H-atoms which were not removed by the 
probe increased from 1 to 5 % as the linear flow velocity was Increased from 
500 to 1700 cm sec -1 . Since kinetic data were collected mainly at velocities 
near 500 cm sec -1 , the probe is thus about 99 % efficient in recombining H-atoms. 

B. Energy Accommodation 

The energy accommodation at the detector, i.e. the question whether 
the difference in electrical power between atoms off and on as measured by 
the Kelvin bridge can be set equal to the power due to the full recombination 
of the H-atom flow, was tested by comparison with H + N0C£ titrations. 23 
Here the quantitative, fast reaction H + NOCfc HC1 + NO removes H and makes 
NO, so that the intensity of the red HNO* emission first rises and then 
decreases to zero with increasing NOCfc additions. The concentration of NOCH 
added at the end point, when the glow is barely extinguished, is exactly 
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equal to the H-atom concentration if the titration reaction is considered 
infinitely fast. As long as H is in excess, I = I 0 [H][NO] = I 0 [H][NOC£] added 
and [H] + [NOC£.] added = [H] 0 , the original H-atom concentration to be 
determined. Rearranging, we obtain T ~ Thnrfl = ^H] 0 “ [NOCi ] a dded and 

a plot of vs. [NOCsJadded should be linear and have an x-axis 

intercept at which [NOC£] added equals [H] 0 as is shown in Big. 6. These 
titration plots have large probable errors for several reasons such as the 
difficulty of controlling and measuring pressures and flows of highly 

T 

corrosive NOCji and extrapolation to zero ~ For each of the 10 

experiments listed in Thble 1, [H] was measured independently by the catalytic 
probe method and by NOCK, titration , with the result that [Hlprobe = (° -99 ± 

0.08) [H] 0 , showing that energy accommodation at the probe is quantitative 
and ruling out the possibility that vibrationally excited H 2 may be formed 
at the catalytic surface and carried downstream unrelaxed. 

C. Flow Characteristics 

Calculation of kinetic rate constants from measurements of absolute 
atom concentrations as a function of distance along the flow tube requires 
an understanding of the flow parameters , of the perturbation introduced by 
the probe, and of all other simultaneous atom loss processes, particularly of 
the rate of wall recombination. For the first of these three problems, the 
characterization of the flow, one needs to justify the use of the simple 
one-dimensional plug flow approximation by showing that the viscous pressure 
drop, radial, and axial diffusion effects are all negligible. Under typical 
conditions such as pressure ~ 6 torr,.v ~ 500 cm/sec, the measured (and 
calculated ) press ure drop was < 1 %. Radial concentration gradients were 
•estimated by the relation Ac/c - r 0 2 (k^ + 3k w )/8D, where k v and are effective 
first-order constants for volume and wall removal of atoms and D is the diffusion 
coefficient. 24 Typically, Ac/c was calculated to be about 1 % showing this 
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effect to be negligible, too. Axial diffusion leads in a one-dimensional 
approximation to an underestimate of the true, effective first-order rate 
constant such that k 1 = k* bs (1 + k* bg D/v 2 ), where v is the average flow 
velocity in the tube. If the quantity k^g D/v 2 <-c l > the effect is 
negligible. Under typical conditions, k^ bg <. 10 sec -1 , D - 100 to 200 cm 2 
sec -1 , v ~ 500 cm/sec, the above quantity is <, 0.01. The assumption of sixtple, 

j j 

one-dimensional plug flow, = v is therefore justified. 

The second problem, perturbation by the probe, can be estimated by 
the approximate expression 21 

(i . slj = . (72 + i! k 1 D )' /2 ^ 

where c* is the actual concentration x cm upstream of the probe and c is the 
hypothetical unperturbed concentration in the absence of the probe. Substi- 
tution of experimental parameters shows that for xi2 cm the perturbation 
was < 1 % so that probe effects could also be neglected. 

D. Surface Recombination 

Through proper selection of system parameters and by operating under 
conditions of negligible pressure drop, axial diffusion, and radial concen- 
tration gradients, the kinetic rate expression is reduced to 

- v » k, [B] + 2 kfjBO* [M] , 

i.e. the sum of the first order wall loss and the homogeneous tenn. Here 
we assume that recombination with M = H can be neglected, since under most 
of our experimental conditions [H] <_ 1 % and recent estimates 12 indicate 
that k R H A R Hz <L 2 near 300°K. 

The relative importance of these remaining terms has been the subject 
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of some debate, and failure to include correction for possibly large wall 
losses has led to serious error in reported rate constants. In our system, 
the surface's of the downstream dissociation and the cylindrical flow tube 
were poisoned with syrupy phosphoric acid and pumped until the residual flow 
of volatile species due to evaporation was less than 10“ 7 of the total 
flow under typical experimental conditions. Phosphoric acid has often been 
used for reducing surface loss and the values we obtained for the surface 
recombination coefficient were measured repeatedly and found to be between 
10 -5 and 10 -6 . In order to measure the first order wall loss, one should 
Operate under conditions where the wall loss tern is greater than its homo- 
geneous counterpart, i.e. where 

'“wall > 2k M M CM] . 

At 0.3 torn and [H] - 8 x I'O 13 atoms/cc, the pseudo first order term due 
to homogeneous recorrtoination is reduced to about 0.01 sec" 1 and the loss 
due to heterogeneous reaction can be measured directly. However, in these 
experiments , the above mentioned corrections for axial diffusion and probe 
perturbation must be applied. Experiments performed over a period of 
several years showed the wall loss to be small and easily controllable. 

Only infrequent recoating was required and the catalytic efficiency of the 
newly coated wall could be reduced through continued exposure to H/H 2 flows. 

The following is a chronological room temperature evaluation of for a 
freshly poisoned surface measured over a period of 8 days: = 0.81, 0.80, 

O. 67 , 0.42, 0.43j 0.39, and 0.20 sec -1 . The value then leveled off and 
typically was found to be about 0.15 sec -1 corresponding to a surface recom- 
bination coefficient of 1.7 x 10“ 6 . 

Under experimental conditions established 'to measure the homogeneous 
reaction rate constant , i.e. pressures greater thari 5 torr, [H] > 9 x 10 14 atcue/cc. 
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the heterogeneous loss contributed typically 3 % to the observed falloff and 
the values for were adjusted accordingly. Prior to our preliminary pub- 

lication of some of the present results, 25 k w had been estimated less 
accurately from experiments at several torr pressure, and assumed to be 
totally negligible. The small differences between kj^ 2 values reported here 
and earlier 25 is entirely due to the subsequent subtraction of the small 
k^ term. 

The behavior of k^ and probe efficiency at tenperatures below 298°K 
is of considerable - importance. As the NOCJ, titrations could not be carried 
out at 77°K, another method had to be used to check on the consistency of 
probe detection efficiency at low tenperatures. Under conditions of constant 
H 2 and H-atom flow and of negligible reaction loss due to either homogeneous 
or heterogeneous recoirbination, the atom flux was measured at 298°K and the 
tube was then cooled to 77°K. No change was observed indicating constant 
detection efficiency. Also, experimental measurements of at lower 
temperature showed it to be as small as or perhaps slightly smaller than 
its value at room temperature. Since kp more than doubles as T decreases 
from 298 to 77°K, the correction due to heterogeneous reaction decreases 
to about 2 % at 77°K. 

E. Typical Kinetic Experiment 

As in earlier studies involving catalytic probes, 21 a voltage signal, 
simply related to the electrical power required to keep the probe at its 
preselected resistance (i.e. temperature ) , is obtained as a function of 
distance; and two such voltage vs. x traverses, one with atcms off and the 
other one with atoms on, provide an [H] vs. x (and thereby vs. t) scan. Our 
mechanical drive technique, coupled with the self-balancing Kelvin bridge, 
allows us to obtain a complete , continuous [H] vs . x scan in a few minutes , 
whereas earlier Investigators positioned their detector point by point along 
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the flow tube and had to wait at each position for full thermal equilibration 
of the heated detector with its inmediate surroundings.. _ In the present 
arrangement, local heating effects cancel, because the probe, which is~~"-~~ 
automatically held at constant temperature (normally 20 to 50° above ambient), 
is moved uniformly along the tube and thereby represents a local perturbation 
of the tube which is constant in time. Ibis cancellation of potential errors 
was repeatedly confirmed by showing that [H] vs. x plots were Independent of 
probe speed. 

A typical set of x-y recorder traces, one each with atoms off and on, 
is shown in Fig. 7 and reveals the following interesting features: A rise 

in signal at the downstream position (near x = 0) due to local heating by 
the stationary probe near its starting position, because space limitations 
made it inpossible to retract the probe farther downstream. As the probe 
moves upstream, away from its slightly heated initial position, the signal 
rises because more electrical power is required to keep it at its preselected 
resistance (temperature) ; a "dip” near 50 cm, which was entirely reproducible, 
corresponded to a rotation of the probe due to twisting of the long lead 
wires at the extreme downstream end of the flow tube extension. As it was 
present equally in both the atoms off and on runs at the same tube position , 
it had no effect on [H] vs. x plots. 

Making use of the results of Sections III A, B, and C, we obtain 
the following expression for [H] from the measured difference in electrical 
power between atoms off and atoms on at any point along the flow tube : 

[H] (atone cnr 3 ) * P .. (tprr) ft°° ] 2 V off ~ V - on 

4 .184 x 760 x 4H x F l 100 + 

l 6eC 

where AH = 52.1 kcal/gatom H, Rg (ohms) is the preselected resistance, normally 
near 2.5 ohms, and the V*s (volts) are the voltages across the 1 ohm resistor 
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as read off plots such as Pig. 7- The results of Section III D further 
suggest that data analysis can proceed by simple 1/[H] vs. x plots, and that 
the corrected specific rate; constant, k^, can be calculated from the . \ 

uncorrected k mc = as obtained in Pig. 8, by subtracting a small \ 

f ky 

pseudo-third-order wall correction tern, k^ = g [m ] [ h ] — » where [H] av is 
the average value of [H] for that experiment. Although this last correction 
tern artificially and wrongly expresses surface recorrbination as an equivalent 
third-order reaction, such a procedure is justified because k^ is never 
larger than about 5 t of ^ and the rigorous computer analysis of such 
data using simultaneous first and second-order (third-order with [M] constant) 
terns would be unnecessarily awkward and tine-consuming. Typically, for 
experiments in pure Hg at 298°K and 6 torr with \% H and a factor of two 
decrease of [H] along the tube, [H] av ~ 1.5 x 10 15 cmr 3 , and with ky - 0.15 sec -1 , 
ky - 0.25 x 10 -33 cm 6 sec -1 or about 3% of measured ^ ~ 8.3 x 10 -33 cm 6 sec -1 . 

IV, Results and Discussion ■ 

A. k R Hz at 296°K; Precision and Accuracy 

The measurement of the H + H + rate constant at 296°K was first 
carried out over a wide range of pressures (~ 2 to 15 torr), %H. (0.2 to . 7 ) 
and flow velocities (200 to 1300 cm/sec) to obtain a consistent and reliable 
value and to establish the kinetics of the process as cleanly second-order 
in [H] and first-order in [Hg]. Although [H] decreased by only a factor of 
2 to 3 in a. typical experiment (Fig. 8), which is insufficient to prove 
second-order behavior, the constancy of the second-order rate constant at 
constant [Hg] for widely different [H] 0 , together with precisely linear 1/[H] vs. 
x plots proved the [H] 2 dependence of the rate. The linear [Hg] dependence 
of the rate is shown in Fig. 9 in which the corrected second-order rate constants 
of 37 experiments at 2 to 15 torr are plotted vs. [Hg]. The slope of the 
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straight line plot is 8.0 ±. (o = 0.2) x 10" 3 3 cm 6 sec -1 . In excellent agree- 
ment with a value of 8.1 ±0.4 x 10“ 3 3 for 73 experiments at 5 to 7 torr at 
0.2 to 1 . 1 % H. The continual verification of this rate constant was then 
used as a secondary standard in experiments at different temperatures , with 
different M-gases, and with D 2 in place of H 2 , over a period of three years 
during which the probe was replaced once, the dissocldtor several times, 
the tube was re-coated several times, and flowmeters and transducers were 
replaced and recalibrated. Serious deviations from 'the above value of ,k R ^ 2 
invariably led to the discovery of an experimental malfunction such as a 
leak or incomplete atom recombination on the probe (which was heated .to 
dull redness in an H-atom stream at regular intervals)-. In all, about 200 
such experiments were performed over a period of three years . They can be 
grouped as follows: - 100 runs (early work, T dependence, some D 2 work) — 

8.1 + 0.4 x 10~ 33 -cm 6 sec -1 ; 25 runs (mainly M = He work) - 8.1 ‘± 0.3; 

31 runs (mainly M = Ar work) - 8.3 + 0.3; 22 runs (new 1 detector) - 8.1 *+ 0.5s 
An overall value of 8.1+ 0.4 x 10 -33 cm 6 sec -1 seens fully established Wherf* 
the 5 % uncertainty equals a single standard deviation from the mean. It is 
difficult to estimate the contribution of systematic errors arising ffrom 
the absolute accuracy of pressure and flowrate measurements, of flow 'tube 
dimensions ,' and of electrical components and meters. If these add an uncer- 
tainty of about 8$, we may estimate a total standard deviation of + 1055. 

As a f lowtube atom reaction rate measurement it appears to be one of 'the 
most precise and accurate. 

: ’B. Results for k R Hz at T < 298°K and for k R He and at 77° and 298°K. 

Experimental results at 10 temperatures below 298°K are summarized 
in Table 2. Each set of low tenperature data took about one day to collect 
and was nested between check runs at 298°K. Temperature equilibration at 77°K 
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was far more easily achieved than at intermediate temperatures where it required 

a stable balance between liquid nitrogen cooling rates and thermal losses. 

In achieving the flat, constant T vs. x profiles in the flow tube which are 

a prerequisite to a successful kinetic run, the catalytic probe was used < 

as a sensitive, scanning resistance thermometer (with the atom source off) 

and T profiles were recorded (at >. 5 minute intervals to dissipate the effect 

of probe heating) until they were sufficiently constant in both x and t. 

The large number of experiments at 77°K reflects the ease of attaining that 

temperature as well as the need of collecting k^ 2 data at 77°K whenever 
He Ar 

kpj or k R data were collected at that temperature. 

Experiments with M = He or Ar performed with those gases added either 
upstream or downstream of the dissociator gave identical results. In any 
series of runs, data were collected for pure H2 and for anywhere from two to 
six different mole fractions of He or Ar, %, between 0 and 0-95- On the 
reasonable assumption that the measure rate constant, k T = X^ k Hz + k^, 
the extrapolated linear k T vs. Xj^ plot gives k M at X^ = 1. Such a plot for 
M ® He at 298°K is shown in Fig. 10. Assuming that errors may be slightly 
correlated in time, each series of X^ runs was analyzed to obtain a k^Vk^ 2 
ratio and these ratios were averaged as shown in the first two columns of 
Table 3. The ratios were then multiplied by the best, averaged k^ 2 at that 
tenperature to give k^ as shown in the last two colurms. k He is seen to be 
13 ± 6? smaller than k^ 2 at 298° and 35 ± 9 % smaller at 77°K whereas k^ 1 * is 
1^ + 9? larger than k^ 2 at 298° and ^8 + 30? larger at 77°K which indicates 
an increasing temperature dependence from He to H 2 to Ar. 

C. Deuterium Atom Reconbination at 77° and 298°K. 

Measurements of the D + D + D 2 rate constant were stralght-fotvrard 
and entirely analogous to their H counterparts , but the hi#i cost of D 2 made 
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It necessary to use H 2 for all Initial conditioning including check runs of 
k H2 and k w , then to switch to D 2 at low flowrates for about 20 minutes, 
then to measure ky for D-atoms at low pressure , and finally to collect data 
for the homogeneous recorrbinatiDn at 4 to 10 torr. k^ was found to be about 
equal to its value for H-atoms, but since k^ 2 is 25$ smaller than k^ 2 , 
the surface correction term was 5 % , slightly larger than for its H 2 counter- 
part. More than 30 experiments at 298°K gave k^ 2 = 6.1 £0.3 x 10“ 3 3 cm 6 sec -1 
or a k D2 ratio of 0.75 £0.05. The result of 8 experiments at 77°K was 
15.1 £ 1.0 x 10“ 33 cm 6 sec -1 or a \P 2 A^ 2 ratio of 0.82 £ 0.11, 

D. Comparison with Earlier Experimental Results. 

Although much of the early work, especially that of the 1930's, 
was carefully done, experimental conditions were such that it is inpossible 
to unravel the contributions to the observed kinetics due to surface 
recombination, thermal effects, and M-effects (particularly for large 
[H]/[H 2 ]). This leaves only the work of Larkin and Thrush 10 * 11 and that of 
Bennett and Blackmore 12-14 for comparison. For M = H2, Larkin, and Thrush 
reported 9-4 ± 1.4 x 1CT 33 cm 6 /sec at 293°K (9-3 + 1.4 x 10 -33 at 298°K) 
which is in satisfactory agreement with our 8.1 £ 0.4 x 10 -33 . Bennett 
and Blackmore first reported 12 3*7 ±0.5 x 10“ 3 3 in a clean quartz tube at 
2 to 9 Torr pressure, then 4.7 ± 0.7 x 10 -33 at 50 to 250 torr in the quartz 
tube 13 , and finally 5-8 + 0.7 x 10 -33 at 1.4 to 5.9 torr in a Teflon-coated 
flow tube. 14 Surface recombination was taken into account in all three 
studies. Moreover, these authors found a strong negative correlation of kp 
with k w , an effect never encountered by us. The upward trend of kp 
amounts to 57$ going from ref. 12 to 14, yet is termed 'good agreement', but 
seems puzzling as do two other observations: (1) the unusually large degree 

of dissociation of highly purified H 2 in microwave discharges reported in 
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ref. 12 which runs counter to the experience of several other groups and 
s ugg ests ; that the electron spin resonance measurements substantially over- 
estimated? [H] in these experiments and that k R was therefore too low; 

(2)i the k R D > D 2/k^ ,H 2 ratio of 1.3, in substantial disagreement with our 
0.75 and /with theory . 26>27 

Poor agreement also abounds for' M-effect at 298°K. For k R ^ e /k R **2 
Bennett and Blackmore 13 report 1.3 compared to our 0.87. For kj^Vkj^ 
they find 1.7, Larkin and Thrush 10 0.66, and we 1.14. The Bennett and 
Blackmore M-data 13 refer to their high pressure flow tube experiments 
whose interpretation in one dimensional plug flow .using the average flow 
velocity and neglecting the parabolic velocity distribution as well as the 
slow radial diffusion probably' introduces large errors as Poirier and 
Carr's 28 analysis suggests. 

Only Larkin and Thrush 10 report some teirperature dependence results 
near 300°K. Expressed in T-n form, their data at 213, 293, and 3^9°K 
give n = 0.68, in satisfactory agreement with n = 0.8l based on our kj^ 
values at 77 and 298°. A similar, empirical fit of our two kj^ e points 
gives n = 0.40 and a weighted least squares fit to all of our k R H 2 data of 
Table 2, except for the I =,77°K point which falls below the good straight 
line on a log-log plot, . gi-ves n' : = 0.68. Two point fits of k R 2 using only 
77 and ;298°!.:and of k R ^2 give 0.6l and 0.67, respectively. 

E. Comparison /with Theory 

; It i has been the multiple curse laid on our detailed understanding 
of ' atom recorrbination processes in general, and that of H in particular, 
that a ' variety of simple- minded models of the process give a very reasonable 
account of the magnitude of the rate constant in zero order approximation, 
but that very little improvement is achieved when the models are refined 
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and are then expected to describe the process In detail. Moreover, the 
experiments are difficult and the total range of the quantity, kp, whose 
dependence on T and on M forms the principal criterion on which the 
models are evaluated. Is very small, indeed. Thus, in the present work 
a factor of four temperature change produces only 1.7 to 3 fold changes in' 
Icr, and M-effect ratios range from a minimum of 0.65 to a maximum of 1.48, 
whereas a single standard deviation of a given k^ may be 5 to \5% and that 
of ratios 7 to 20 % . Considering these limitations it is surprising how 
much indirect information can be extracted from the present data. 

The principal questions are as follows: A. Is it possible to 

determine if and to what extent a given H + H + M recombination proceeds 
via the energy transfer (ET) mechanism, the radical-molecule -complex (RMC) 
mechanism, or both? B. For the H-atom recombination at low temperature 
(<_ 300°K) can it be shown whether the process proceeds principally via 
orbiting resonances of H 2 (and HM), i.e. by quantum mechanical tunneling 
through the rotational barrier? If these questions can be answered, one 
may then hope to extract from the data some information on the, interatomic 
potentials . 

Benson and Fueno 29 calculated kp^ r in the ET approximation using a 
steady-state step ladder model with either Morse or Iennard-Jones functions 
characterizing the vibrational states near the dissociation continuum. In 
magnitude, their Lennard-Jones results of 1.00 or 1.05 x 10~ 32 cm 6 sec -1 
(depending on whether gas kinetic or Sutherland deactivation cross sections 
are used) are much closer to the experimental 0.92 x 10” 3 2 . than their Morse 
results of 5.0 or 3-7 x 10 -32 , but neither model comes close to giving the 
correct temperature dependence as they give n < 0.2. 

Keck’s classical variational phase space theory was recently applied 
to H-recombination by Shul and Appleton 20 (SA). The orbiting resonance theory 
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was proposed by Roberts, Bernstein, and Curtiss 18 (RBC), recently modified 
by Whitlock, Muckerman, and Roberts 30 (WMR), and extended to include con- 
tributions by the RMC mechanism by Pack, Snow, and Smith 31 (PSS). The 
comparison of our data with these theories is shown in Figs. 11 (I^H^) , 

12 (M=He), and 13 (M=Ar) . It is best, in posing questions A and B to 
examine all three figures simultaneously. The original orbiting resonance 
theory 18 (RBC) predicts the correct magnitude of k R ^ to better than a factor 
of two, but shows much too small a temperature dependence due to its following 
shortcomings: (1) Too crude a model for the calculation of rotational 

relaxation cross sections, e.g. assuming metastable H 2 * to be a rigid rotor; 
(2) Neglect of the RMC mechanism which could lead to underestimates Of kp^ 
and of all k R M at the lowest temperature ; (3) Neglect of metastable, non- 
resonant states of H 2 and of free states with energy above the rotational 
barrier; (4) Neglect of departures from equilibrium. Correction of (1) 
above was made by WMR for M = H 2 using classical trajectory calculations. 

This produced improvement but still predicted too flat a T-dependence 
at low temperature. A second calculation by WMR, not shown in Fig. 11, 
followed SA and PSS in using an attractive H - H 2 potential with a 38°K 
well depth (whereas the WMR curve in Fig. 1 had used a repulsive H - H 2 
potential) to calculate the Hg - H 2 relaxation cross section. It duplicated 
the shape of the experimental curve, but is everywhere too high by about 
20 to 40%. The classical SA calculation also reproduces the shape of the 
k R ^ vs. T curve, but falls 20 to 30 % below the experimental points even 
though it contains a substantial contribution from the RMC mechanism and 
uses an H - H 2 potential which is probably too strongly attractive. If the 
potential proposed by Dalgarno, Henry, and Roberts 32 , which supports no 
bound levels (e/k = 15°K), is correct, much of the RMC contribution would 
have to be removed and both the magnitude of the SA rate constant and its 
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tenperature dependence would then be Incorrect. Low energy (~20°K) beam 
scattering experiments on H + H 2 would be particularly welcome, since 
data which are presently becoming available 33 are in too high an energy 
range to probe the depth and shape of the potential well . The PSS calcula- 
tion extends the RBC theory by including the contribution which M-H 
resonances make to the recombination, i.e. by including the quantum 
mechanical RMC mechanism, and by evaluating the Inelastic collision cross 
sections in the WKB approximation. For M = H 2 , these calculations over- 
estimate kR somewhat near 300°K and underestimate it at 77°K even though 
they use the deep (38°K) attractive well for H-^* It is interesting to 
note that the addition of the RMC contribution (as indicated separately in 
the PSS paper but not in Fig. 11) to the WMR rate constant would produce 
good agreement with our experimental values over the entire tenperature 
range (but less so if the H-H 2 potential has a shallower well). 

For M = He, all theories prodict a flatter kp vs. T plot than the 
admittedly slim experimental evidence indicates. SA's phase space calcula- 
tions are able to reproduce such a steeper dependence (not shown in Fig. 12) 
by using a well depth of 39°K, r e = 3-7ft, and e = 1.87ft -1 in the Morse 
potential for H-He. This probably is much too strong an interaction. 

Recent theoretical 34 and experimental 35 studies indicate a well depth in the 
5 to 10°K range, in line with the now accurately measured He-He well depth 
of 11.0 + 0.2°K 36 . Such a potential does not support a bound H-He state and 
removes the possibility of a RMC recombination pathway as shown by the 
flat maximum in k R at 100 to 150°K calculated by PSS. The older RBC calcula- 
tion comes closest to fitting the data, but in view of its approximations, 
this is probably not significant. * 

For M = Ar, all recent calculations fall below the experimental 
points and only the RBC curve has a maximum. PSS report a smaller IMC contri- 
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bution In the M = Ar than in the M = H 2 case which is surprising, particu- 
larly since a sizable well depth of ^7-7°K is used for H - Ar (r e = 3*36^) • 
For a reasonable choice of potential, the classical SA curve has the 
correct shape , but again falls well below the experimental values . 

The k R D A R ^ ratio is in good agreement with both classical and 
quantum calculations. The latter 27 also predict an increase of this ratio 
with decreasing temperature which is marginally supported by our observed 
rise from 0.75 ± 0.05 to 0.8l + 0.11. 

Tb return to the questions posed earlier: it seems that the M = Ar 

data and the M = H 2 results at the lowest temperatures are not easily 
rationalized without appeal to the RMC mechanism, but that it plays a minor 
role under all conditions investigated here except, perhaps, for M = Ar 
near 77°K. It would be desirable, therefore, to extend recombination rate 
measurements to still lower temperatures , particularly for M = He . 

On the question of quantum vs. classical models a clear decision 
can not be made. The resonance theory Is physically more reasonable and 
has shown good agreement with experiment at an earlier stage in its develop- 
ment, i.e. with less refinement and for more realistic choices of interaction 
potentials. More rigorous calculations should be carried out to ascertain 
the role of non-resonant as well as continuum states as recombination 
intermediates . The present experimental work Is continuing with measurement 
of other M-effects including, hopefully, M = H and of the ortho-para ratio 
of recombined H 2 in an attempt to resolve question B directly 37 . 
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Table 1 


Comparison of [H] Measurement by NOC£ . Titration 
and by Catalytic Probe 


Total Pressure 
(Torr) 

[H] 

( 10 15 atoms/cc): 

. ^probe 
(lO 1 ^ atoms/cc) 

Mprobe/M 

- i%) 

• 2.12 

0.86 

0.86 

93 

1.29 

1.32 

1.17 

89 

2.22 

1.28 

1.19 

93 

1.30 

' 1.29 

1.26 

97 

1.00 

1.27 

1.28 

101 

1.34 

1.25 

1.32 

106 

1.35 

1.59 

1.48 

93 

1.35 

1.3^ 

1.49 

111 

1.94 

1.54 

1.58 

103 

1.95 

1.51 

1.63 

107 



Mean Value 

= 99 ± 8 




Table 2 


kp 2 as a Function of Temperature 


.Temperature 

(°K) 

Number of 
Experiments 

k R x .10 3 3 (cm 6 /sec) 

-uncorrected for 
wall loss 

k R x 10 3 3 (cm 6 /sec) 

corrected for 
wall loss 

298 

> 175 

8.3 + 0.4 

8.1 ± 0.4 

275 

7 

8.5 ± 0.5 

8.3 ± 0.5 

253 

■6 

9.6 + 0.3 

9.4 ± 0.3 

231 

6 

10.0 ± 0.8 

9.8 ± 0.8 

210 

6 

10.8 ± 0.9 

10.5 t 0.9, 

200 

5 

11.3 + 0.2 

11.1 ± 0.2 

180 

6 

11.7 + 0.2 

11.4 + 0.2 

160 

8 

12.2 + 1.0 

12.0 + 1.0 

140 

7 

13.5 ± 0.9 

13.3 ± 0.9 

120 

8 

15.6 + 1.6 

15.4 ± 1.6 

77 

72 

18.9 + 2.2 

18.5 ± 2.2 
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Table 3 

M 

Recombination Rate Constants, kp, at 77 and 298°K 


M 

. . M 

298°K 77°K 

kp x 10 33 
(cm 6 /sec) 

298°k 

77°K 

»2 

1.0 

1.0 

\ 

8.1 ± 0.4 

18.5 ± 2.2 

He 

0.87 ± 0.06 

0.65 ± 0.09 

7.0 ± 0.4 

12.0 ± 1.5 

Ar 

1.14 ± 0.09 

1.48 ± 0.30 

9.2 ± 0.6 

27.4 ± 4.6 

d 2 (D) 

0.75 ±0.05 

0.82 ± 0.11 

6,1 + 0.3 

15.1 + 1.0 



28 


LEGEND FOR FIGURES 

Fig. 1 Diagram of apparatus 

Fig. 2 H-atom source 

Fig. 3 Catalytic probe support 

Fig. • self-balancing Kelvin bridge 

Fig. 5 Probe collection efficiency determination 

Fig. 6 NOCil titration conparison ' 

Fig. 7 Typical recorder traces for measurement' of [H] vs . x 
Fig. 8 1/[H] vs. x plot 

Fig. 9 Plot of second order H-atom decay rate constant, k 1 * vs. pj^ 

Fig. 10 Plot of total third order rate constant, k T vs. 

Fig. 11 k R H 2 vs. T. J - this work; RBC - ref. 18; SA - ref. 20; WMR - ref. 30; 
PSS - ref. 31 

Fig. 12 k R He vs. T. symbols as in Fig. 11 

Fig., 13 kj^ vs. T. symbols as in Fig. 11; A - ref. 10 
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